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We investigate how equilibrium entanglement is manifested in the nonlinear response of an A''- 
qubit system. We show that in the thermodynamic limit the irreducible part of the nth-order 
nonlinear susceptibility indicates that the eigenstates of the system contain entangled (n + l)-qubit 
clusters. This opens the way to a directly observable multiqubit entanglement signature. We show 
that the irreducible part of the static cubic susceptibility of a system of four flux qubits, as a function 
of external parameters, behaves as a global 4-qubit entanglement measure introduced in Ref . [20l| . 

We discuss the possibility of extracting purely-entanglement-generated contribution from the gen- 
eral multipoint correlators in a multiqubit system. 



I. INTRODUCTION 



The "equilibrium" , or "thermodynamic" , entanglement in spin systems is a quantification of the degree of non- 
classicality of the eigenstates of the system (e.g. [l|, 0, S S @)- In quantum phase transitions it can be regarded 
as a specific order parameter [g', '?]. It distinguishes between the classical correlations, produced by interaction, and 
the purely quantum ones. Both types of correlations are contained in multi-point Greens' functions of the system. 
However, there is no simple recipe to separate them. This problem is exacerbated by the absence of a convenient 
measure of entanglement for a mixed state of more than two subsystems. Quantum state tomography can reconstruct 
all qubit states of a general spin system, but the number of required operations grows exponentially with A'^ 0]. 

Various entanglement measures were expressed through spin correlators in Observable "entanglement wit- 

nesses" [ll| or "entanglement estimators" [l^l were proposed as well. In general, a number of measurements is 
required to restore the measure of entanglement in equilibrium, though not as large as in the case of quantum state 
tomography. 

Equilibrium entanglement is particularly important in the context of adiabatic quantum computing (AQC) p^ . 
From the equivalence of AQC and standard quantum computing [13, fTi| and the polynomial equivalence of quantum 
and classical computing in case when there is no global entanglement|l(| we infer that global entanglement of the 
ground state is necessary for an efficient AQC. The equilibrium state of a system is a sum of projectors on its 
eigenstates. It is therefore desirable to be able to check that at least some of these project on globally entangled 
eigenstates. 

The investigation of equilibrium entanglement allows us to shift the focus from the state to the Hamiltonian (e.g. its 
entanglement-generation ability[l7i). This opens a new venue of investigation. Suppose the Hamiltonian of a quantum 
system can be restored without free parameters from a series of measurements (which means that the behaviour of 
the system, including level anticrossings etc, is quantitatively described under the assumption that it is in a thermal 
equilibrium mixed state produced by this Hamiltonian). Then all the eigenstates of the system can be determined, 
and their entanglement measures can be computed. 

The experimental investigation of equilibrium entanglement along these lines was performed for the sytems of 
twop^ and fourflo', '2^ flux qubits. It relied on continuous measurement of the magnetic susceptibility of an A^-qubit 
system. The qubits were controlled locally; the global response was measured, and all the operations were slow on 
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the scale of characteristic qubit times. The quantification of entanglement in both cases was done "post-mortem" and 
showed a high degree of entanglement for the ground and first excited states. 

Remarkably, the results of [la] also contain a directly observable, "in vivo" 2-qubit entanglement signature, so- 
called IMT deficit (a difference between the sum of signals from two qubits driven through their degeneracy points 
separately, and the signal when they pass the co-degeneracy point simultaneously; the abbreviation comes from the 
measurement technique employed). It appears in the linear susceptibility, which is proportional to 

Rp,-Y.^pyi\q){qy-\p), p^q. (1) 

Here \q) are the eigenstates of the Hamiltonian, and (t| is proportional to the magnetic moment of the jth qubit 
in z-direction. (The actual formula (e.g. Eqs.(2,3) of Ref. [la| : see also (2lj) takes into account asymmetries in 
couplings of qubits to the external field source; while necessary for the quantitative match to the experiment, it is 
not conceptually important.) If the eigenstates are not entangled (that is, do not contain components differing by 
two qubit flips in the z-basis), part of this term (with i ^ j) disappears. The only requirement for the observation of 
this signature is, as we stated above, the ability to manipulate qubits separately and measure the total response of 
the system. 

The linear response measurements will not provide a directly observable signature of equilibrium entanglement with 
N > 2. The expressions like ([1]) arise in the lowest-order perturbation theory and therefore contain only first-order 
commutators between the qubit operators (coming from the expression for the observable, and from the perturbation 
Hamiltonian). Therefore Rpq can only catch the absence of 2-tangled components, but is insensitive to the higher- 
order (e.g. GHZ-type) states. This is why the appearance of globally entangled eigenstates (| Titi) i I itit)) in 
could be only inferred from the "post mortem" analysis. 

In order to obtain information about entanglement with N > 2 without reverse-engineering of the Hamiltonian, 
one requires A^-point correlators. In classical case, they can be extracted from the noise fluctuations ("noise of the 
noise"), but it is not clear yet how and whether this can be generalized to the quantum case. A more immediate 
approach is to measure nonlinear response of a multiqubit system. In practice, at least quadratic and qubic effects 
should be observable, thus providing signatures of 3- and 4-qubit equilibrium entanglement. 



II. ENTANGLEMENT AND HIGHER-ORDER CORRELATORS 



First consider an A^-point correlator of single-qubit spin operators, Sk{tk), acting on qubit k at the moment tk, 

CN{{k,tk},\p)) = {p\Sl{h)S2{t2)...SN{tN)\p), k^l....N, (2) 

where \p) is some eigenstate of the Hamiltonian. The sequence of operators in ([2]) generally flips N spins in the state 
\p). Therefore the correlator disappears unless \p) is a superposition of components differing by N spins. In general 
case, one expects that Cn should disappear or at least diminish in the absence of entangled states involving at least 
N qubits. 



A. Equilibrium 7V-qubit entanglement 

In a system of M qubits, every eigenstate of the Hamiltonian can be written as a product of states defined on qubit 
clusters ri^*'' of no more than M qubits each (Figllji), 

in such a way that the qubits from different clusters are not entangled. Here Ai is the set of all qubits in the system; 
\A4\ = M . Since every time we deal with a pure state, this operation does not produce any problems. The complement 
of a cluster we denote by = A4\il. 

For the following, we need to impose on the Hamiltonian the 

Russian-doll condition (see FigH]): If in any eigenstate \p) qubits a,b E rji^"*, while qubit c ^ ili^', then 
there is no such eigenstate \q) and cluster \ that a ^ ^^r \ while 6, c S 
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Equivalently, it can be stated as follows: if any two clusters fli'^'' have a common qubit, then one of these 

clusters must be contained in the other (like Russian dolls): 

n ^ ^ (n'-p^ c v (ni'^^ c n'f ')^ . (4) 

The Russian-doll condition allows us to uniquely partition the M-qubit system in equilibrium into mutually non- 
entangled clusters (maximal clusters) flj, 



^ — n 



(5) 



each containing no more than M qubits (see Fig|TJ:,d). This is a physically plausible picture, and we conjecture 
that most iV-qubit Hamiltonians satisfy the Russian-doll condition. Establishing the limits of applicability of this 
conjecture will be the subject of further research. 
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FIG. 1: (Colour online.) Russian-doll condition and partitioning the system into maximal clusters, (a) An eigenstate \p) of an 
M-qubit system can be presented as a product of functions, which depend only on states of qubits belonging to the clusters 



In other words, an eigenstate determines the partition of the system into mutually non-entangled clusters, (b) 



Two sets of qubit clusters, for eigenstates \p) and \q), are shown. The Russian-doll condition is not satisfied, since qubits a and 
h are entangled in \p), but not entangled in \q), while qubits h and c are not entangled in but entangled in \q). (c) Russian- 
doll condition is satisfied for three eigenstates. The entangled clusters for them are marked by solid, dotted and dashed lines 
respectively. The largest "Russian dolls" are the maximal clusters (clusters containing the largest number of qubits, mutually 
entangled in some eigenstate of the system), (d) The unique partition of the system of Fig.(c) into maximal clusters (clusters 



In a system satisfying the Russian-doll condition, the existence of a fc-qubit maximal cluster is equivalent to any 
state of the system being fc-producible (see Refl2^ We can therefore say that equilibrium A^-qubit entanglement 
exists in a system oi M > N qubits, if at least one maximal cluster contains at least N qubits. Indeed, then in 
equilibrium the density matrix of the system is a sum of projectors on eigenstates, at least one of which supports 
iV-qubit entangled clusters. 

Our goal is to find an observable signature of equilibrium A^-qubit entanglement, that is, an observable physical 
quantity which will be zero unless the system contains maximal clusters of size N . Let us check the handwaving 
arguments we made around Eq.([2]) and see whether certain A^-point correlators can be used this way. 
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B. Analysis of irreducible correlators of one-qubit operators 

Returning to ([2|), one can easily express it as a sum over the internal labels P2,P3, ■ ■ -Pn of the correlators (pi = p) 
CN{{k, \pk)}) = {pi\Si\p2){p2\S2\P3) . . . {pn\Sn\pi), k^l....N, (6) 

with the appropriate energy denominators, similar to ([1]). (This is done by the repeated insertion in ([2]) of the closure 
relation, i = \p){p\-) 

Now let us concentrate on the irreducible correlators, that is, the ones where no two states \pi), \Pj), i ^ j coincide. 
The term reflects the following property of the correlator ([6]). Obviously from the definition, CAr({fc, \pk)}) is invariant 
with respect to the cyclic transpositions of the labels pk- If some of these labels appear more than once, then 
C7v({fc, \Pk)}) is also invariant with respect to shorter cyclic transpositions in substrings, which end in the identical 

PkS. 

For an irreducible correlator c'^''({fc,p}) there are two possibilities: either all N qubits, on which act the operators 
Sj, belong to the same maximal cluster 51, or they belong to different maximal clusters. 

1. Same maximal cluster 

In this case the states of the system can be written as \n) = |n)o ® |n)Q etc, where Cl is the complement of 51: 

Vtyj Cl ^ M. 

(This follows from the definition of entanglement and the fact that in any eigenstate of the Hamiltonian 17 and 57 do 
not contain mutually entangled qubits.) We therefore obtain 

<({fc,ri) = (• • ■)5j,,j,M5p.pM ■ --^.p^m (7) 

(the Kronecker symbol Spq{ft) means that the states p and q coincide on the set 51). Eq.([7]) reflects the fact that the 
states of the qubits outside il are not changed by the operators acting on qubits inside 51, and can therefore only 
differ on 57. The number of states which only differ on a cluster of size |51| = Nfi is 2^*^. Therefore there can be no 
irreducible correlators with N > 2^" (otherwise the indices must repeat). Of course, this is not a serious limitation. 

2. Different maximal clusters 

Now the string of indices 

which labels the qubit operators, consists of 1 < i? < substrings, S — siS2 ■ ■ ■ Si?. In each substring the indices 
belong to the qubits in the same cluster, and there are R such clusters involved: fli, 5I2, . . . 51^ (see Figl^^). Some of 
these clusters may coincide (i.e. called more than once), with the only trivial restriction: Clusters corresponding to 
adjacent index substrings, e.g. 57^ and 57r±i, must be different; otherwise these substrings concatenate. 

The rth substring (the indices in which run from j[ to jj^ , and j]^ + 1 = j]^^) yields the product of the Kronecker 
symbols, ensuring that on 57^ (the complement of 57^) 

\Pj{) = \Pji+i) = ■■■ = \Pr;,^+i) = bj[+i>- (8) 

To shorten the presentation, we introduce the following definitions: 
Definition 1: Substrings are mutual, if they contain qubit indices from the same maximal cluster. 
Definition 2: A unique substring is the substring, which has no mutuals in S. 

Definition 3: A joint correlator is the correlator, the index string of which contains at least one unique substring. 
The index string of a disjoint correlator does not contain unique substrings. 
Now consider a joint correlator with the unique substring (Figl2)D,c): 



S = S1S2 ■ . . Su-lSuSu+l . . . SR-iSr. 
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FIG. 2: (Colour online.) Joint and disjoint irreducible correlators, (a) The set of all qubits, A4, and maximal clusters called 
by the index string S. (b) The complement A4\Q,u = Qu, oi the cluster Qu- (c) If ilu is unique, then the complement of the 
sum of all other clusters called by the string S, M\{Qs UQpU^lt), contains (d) Therefore the sum of (b) and (c) is the set 
of all qubits, Ai. (e) If Qu is not unique, the rest of the called clusters contains and its complement does not. (f) Then 
the sum of the complements, (b) and (e), is not M, but only Clu (f). 



Then, from Eq.® we obtain: 

(• ■ ■) ^pipj^ (Hi n • ■ • n Qu-i) Sp.^p [flu) Sp ^^ipi {fiu+i n ii„+2 • • • n Ciu) = 
1 1 jj 

(• ■ ■) Sp^^p.^+i {^u u (f^i • • ■ n Ciu-i n n.u+i • ■ • n Qr)) 



(9) 



Since s„ is unique, that is, flu H = 0, the set flu U {fli ■ • • fl riu_i fl flu+i • ■ ■ H 17^) = In other 

words, the states and must coincide everywhere; therefore they are identical: \pju+i) = \pj^). This means 

that a joint correlator cannot be irreducible. 
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Nevertheless, for a disjoint correlator this does not hold (Fig[2F): 

17„ u (f^i • • ■ n Cir) = n^ = M\n,^ c M, 

and we cannot conclude that = The only exception is the case when si and sr, the substrings at 

the beginning and the end of 5, are mutual, and don't have other mutuals. Then, from the cyclic invariance of the 
correlator, these substrings concatenate into a unique substring. 

So far we established that an irreducible correlator of the A^th order does not exist, unless (a) the corresponding 
maximal cluster n contains at least Nq — log2 N qubits, or (b) it is disjoint (its string of indices does not contain a 
unique substring). This falls short of our goal to find a signature of equilibrium A^-qubit entanglement. Nevertheless 
the situation changes when we turn to the global response of the system. 



III. GLOBAL RESPONSE OF THE SYSTEM IN THE THERMODYNAMIC LIMIT 

The global response of the system is determined by the reaction of its total magnetization to the external pertur- 
bation. It is defined through the average operators, i.e. operators of the type 

M 

5 = M-i^5„ (10) 
and contains in particular the irreducible correlators of average operators: 

M 

Cn]piP2...pn ^ ^^^^ X! c'^'{ki,k2,---,kN;pi,P2,---,PN)- (11) 

ki,k2,. ..,kj^ — l 

The latter are expressed through the irreducible correlators we discussed in the previous Section. 

We will show that in the thermodynamic limit, when both the number of qubits M and the number of maximal 
clusters Q — > oo, while the order of entanglement N stays finite, the irreducible global correlator C]^"^ (|lip is zero, 
unless a finite ratio of maximal clusters in the system contain > N qubits each. This makes C]^' a signature of A^-qubit 
equilibrium entanglement. 

In the previous section we saw that there are two types of irreducible correlators, which do not disappear, if the 
system does not contain A^-qubit maximal clusters. We will deal with them separately. Without loss of generality, 
assume that M > N and Q > N. 



A. Contributions to from correlators with repeated indices from the same maximal cluster 

The number of all possible combinations of A^ indices, each running from 1 to M, is Ai^ . The number of combi- 
nations of different indices is M{M — l)...(Af — A^ -|- 1), and the number of combinations with repeating indices is 
— M{M — 1)...(M — A^ + 1) = 0{M^^^). Therefore in the expression (jlip the contribution from the correlators 
with repeating indices scales as 0{l/M). 



B. Disjoint correlators 



There are ^ ^ ^ ^ ways to split a string of A'^ indices in R nonzero substrings. For each substring there 

are Q maximal clusters we can choose from, and there are on average (M/Q)^/^ possible combinations of in- 
dices within a cluster. Therefore the total number of different combinations of indices, counted this way, is 

{M IQ)^ X]^=i ( ^ \ \ ■ other hand, the number of combinations corresponding to joint correlators is 



no less than [M /Q)^ Tlfii^i ^ | J ^{Q ^ {R^ !))■ Therefore the number of disjoint combinations does not 
exceed 
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(12) 



and the disjoint contribution to the average irreducible correlator (fTTj) asymptotically disappears at least as fast as 
0(1/Q). 

IV. NONLINEAR SUSCEPTIBILITY AS A SIGNATURE OF ENTANGLEMENT 

The correlators like (|lip enter the expressions for (nonlinear) susceptibilities. Therefore the latter could be used as 
directly observable entanglement signatures. We introduce the equilibrium 7V-th order signature, 

Z^. = (Z;:,)^^Z;^e-^"/^, (13) 

n 

which is the thermal average of the contributions, Z]^, from the sysem's eigenstates \n) with energies En'- 



Pl P2#P1 P1V-1#PJV-2,PJV-1,---P1 



Here /ip^ = (pl/^^k) is the matrix element of the z-component of the system's total magnetic moment. The functions 
Zat, Z"^ are related to the irreducible part of static nonlinear susceptibility of the system. In the Appendix, this is 
explicitly shown for the quadratic susceptibility. 

The largest solid state structure where entangled states were demonstrated so far contains four superconducting 
qubits^lfl^. This is far from the region of applicability of the thermodynamic limit results of the previous section. 
Still, a comparison with the experimental data can indicate whether our approach is usable. 

We will therefore compare the irreducible part of static (w 0) cubic susceptiblity, Y^'^^(aj ) oc Z4, to the measure 
of global entanglement introduced in Ref,'2^ and calculated for the system investigated in Refligi. 

The global entanglement measure for an arbitrary pure state i\j of an n-qubit system is given by[20l| 

/ , X l/(2"--l) 

\1<|S|<S| / 

Here the product is taken over all possible bipartitions of the system, (S', 5); |5| is the number of qubits in the set 
5, and the prime denotes that if |5| = n/2, we still include every variant only once. The quantity rjsO^), < rjs < 1 
characterizes the entanglement of the set S with the rest of the system, S (see Eq.(3) of Ref.2Q,): 



2l'SI 



^75 



2l^^l 



-{l-tr[(tr5p)2]}, (16) 



where p is the density matrix. Obviously rjg = 1 if and only if S and S are separable. The known Meyer- 
Wallach(23| and Scott [25| measures of entanglement are expressed through arithmetic averages of rys's. Unlike them, 
the measure [20j in Eq. (jl5p is a geometric average and therefore is zero if there is even one qubit separable from the 
rest. 

The expression ([T5|) can also be generalized to the case of a mixed state[2^, but we are here content to use the 
definition (jisp and compare the results with the ground-state value of the 4-th order signature for four qubits, 

^ \^ \^ \^ 4 (M^0piM''piP2M^P2P3^^P30) /-.yX 

' . ^ (^fi - ^o)(^f^ - ^o)(^p3 -E^y ^ 

Pl P27^Pl P3#P2,Pl / ^ . V K., 

which is proportional to the irreducible part of the static cubic susceptibility of the system at zero temperature. Since, 
unlike TZ{^), can take negative values, we will take for comparison its absolute value. 
The system in question is described by the pseudospin Hamiltonian 

4 

F = -^ha«+A,a«]+ ^ J.^frW^W, (18) 

i=l l<2<j<4 
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FIG. 3: (Left) Normalized 4-qubit ground-state entanglement signature j22(/i,2, -^'st)!/ rnax/j,^ \Z2{h2,hT)\ (Eq.((I3), and 
(right) global 4-qubit entanglement measure Ti{tpo) (Eq. (|15|| '). as a function of bias currents in the 4-flux-qubit system of 
RefOa. 

where Ai = 147, A2 = 12, A3 = 163, A4 = 165, and J12 = J34 = 163, Jm = J23 = 155, J13 = J24 = -62 (all in 
mK) were determined from the experiment . The operator /i^ oc X]i=i "'z*''; ^^'^ terms in Eq. (|17p can be 

calculated explicitly. 

The dependence of \Z2\ and TZ{ipo) (ground state global entanglement) are shown in Fig[3]as a function of variables 
IbT and Ib2- These variables correspond to the bias currents through different tuning coils in the actual experiment 
and can be explicitly expressed in terms of the biases Ci in Eq. dTH]) (see Refil9i) . FiglJ] represents one possible section 
of the parameter space {ei}i^=i- The centre of this picture is the co-degeneracy point, i.e. here all qubits in pB]) are 
biased to their degeneracy. The maximum value of TZ for four qubits is approximately 0.95 (Refl20l). To make scales 
comparable, we normalized \Z2\ to its maximum value at a given /f,2. 

One can immediately see that the structure and positions of the maxima of both functions practically coincide, 
though the relative heights of the peaks not necessarily so. This is further illustrated in Figdl where the dependences 
on IbT for a given value of Ib2 are presented. Given that the formulae like (|17p are only supposed to hold in the 
thermodynamic limit, such a similarity between the behaviour of H and Z4 is remarkable and indicates that the 
irreducible part of nonlinear susceptibility can provide information about multiqubit entanglement even in modest- 
sized structures available now. 



V. CONCLUSIONS 

In this paper we concentrated on equilibrium entanglement, that is, on the presence of entangled -/V-qubit clusters in 
the eigenstates of the Hamiltonian in a multiqubit system, and its directly observable manifestations. We demonstrated 
that under plausible enough assumptions (Russian-doll condition), and in the thermodynamic limit, the irreducible 
part of static nonlinear susceptibility provides a signature of the existence of multiqubit entanglement in the system. 
We introduced a series of such entanglement signatures. We also considered a four-qubit system, with the choice of 
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parameters corresponding to the experimentally realized situation. We saw that the structures of maxima of the global 
four-qubit entanglement in the ground state, TZ, and of the fourth-order entanglement signature, Z^, as a function 
of control parameters, are almost identical. This raises an expectation that the method we introduced can provide a 
directly observable evidence of multiqubit entanglement even in systems of moderate size. 

There remain open questions. First, it is necessary to investigate the requirements to the Hamiltonians which 
satisfy the Russian-doll condition. Second, one would like to see explicitly the origin of similarity between TZ and -E, 
and investigate whether one can be expressed through the other in some nontrivial limit. Third, it is desirable to 
develop an experimental protocol, which would allow to extract Z of the third or fourth order from the measurements 
and thereby provide a direct evidence of equlibrium multiqubit entanglement. 
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APPENDIX A: IRREDUCIBLE PART OF STATIC QUADRATIC SUSCEPTIBILITY 

We start with the Liouville equation for the density matrix, 

idtp{t) = [Ha + Hi{t),p{t)], (Al) 



where the unperturbed Hamiltonian of M flux qubits and the perturbation are 

M M 

2 



^ M M 

j = l i<j j = l 

The factors Ij answer e.g. for the inhomogeneous external magnetic fleld distribution in the system. The equilibrium 
entanglement we are going to probe is created by the unperturbed Hamiltonian, Hq. (Of course, the following 
discussion is applicable to any system described by the pseudospin Hamiltonian (|A2[) .) 

Iterating Eq. (|Aip , one obtains the standard expansion of the density matrix over the commutators of Hi and po 
(unperturbed density matrix pQ is time-independent and commutes with Hq, but it does not have to be an equilibrium 
density matrix for Ha). The expansion p{t) = p^"^ + p'^^^t) + p'^'^\t) + . . . leads to the corresponding expansion for 

the magnetic moment of the system, p{t) — tr |p(t) ■ 

One can measure different components of the magnetization and apply several components of the perturbation. 
Here we consider the situation closest to the 2- and 4-qubit setup of H, H, [13], where both the external fleld, and 
the measured magnetic moment are along z-axis. For simplicity we also assume that all Vs, are the same. 

The linear term in the expansion of p,z {t) is related to 2-tanglement and was treated in p^ . l2l| . The first new term 

is ^Jif\t), 

^ 1 /■*' r -1 r°° 1"°° 

f,f\t) =J2 ^ dt' dt"h{t')h{t")tT{p^''^ [[af(i),a|(i')] ,'yk{n]}= / dt' / dt"x..{t-t',t-t")h{t')h{t"lA3) 

i j k—1 ^ — oo J — oo ^ J — oo J — oo 

where the quadratic susceptibility Xzz is 

N 

X..{t-t',t~t")^ J2 -^Oit-mt'-t")tr{p'^'^[[a:it),a^{t')],al{t")]}. (A4) 
In Fourier components we get 



ij.k—l 



^ /IS /I^e-(-+-')*x..(c-,c.')M-)M^')- (A5) 
For a monochromatic perturbation, h{ui) — 2TrTld{uj — ujq), the response is p,^^\t) = e~^*"°*7i^Xzz(i^o, t^o)- 
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We can rewrite (IA4[) using the eigenstates of Hq, Ha\n) = -E„|n), n = 1,2 . . . 2^: 

9 A/ 

x..it -t\t- 1") = ^ E E Pn'^o^i n{n\ [[am, <^m] ' ^nn] \^)^ 



(A6) 



ijk n—1 



where all operators cr^(t) are in the interaction representation. Opening the commutators, using the closure relations, 
taking the Fourier transform and doing the usual algebra, we finally obtain: 



(0) 



Tl = l 



n=l 



^ ^ ^n:pqfn\pq{^:^ ) 
p,q=l 



We introduced the weight factors 



Cn-pq - J2^n\<7^\p){p\a'^\q}{q\al\n) = {n\fi'\p){p\fi'\q){q\^i'\n), 

ijk 



and the formfactors 



(A7) 



(A8) 



LU + LU' - {Eq - Ep) + iQ 



X (A9) 



En + Eg ~ 2Ep 



En + Ep — 2Eg 



(w' - (Eg - En) + iO)(w + w' - (Ep - En) + iO) [lo' - {En - Ep) + iO){uj +L0'- (En - Eg) + lO) 



Obviously, Cn-^pg = 



z ,,z ,,z 



n — ^npq — ^pqn — ^qnp, and fji-jmi.^ 7 ^ ) 0- 



.pq — r T V Xlr. 9" ~ "^ripq — '^pqn — '-^qnp 



In the setup of [18. ,1191, 120| the frequency of the LC-tank, at which the susceptibility is measured, is orders of 
magnitude less than the tunneling and coupling frequencies 20 MHz vs. ^ 1 GHz). The interlevel spacings are at 
least ~ A/M. Therefore if w 4; A we can use the static limit: 



^0 = lim,E^"P9-/'«;P'3(^'^)' 

a;— >0 ^ — ' 



where means that terms with Ep = Eg ^ En are excluded. Separating the contributions with both Ep ^ En, Eg ^ 
En yields 



EE 



'3C„ 



[Ep — En) {Eg — En) ^ {Ep — En) 



pq 



6EE 

p q<p 



Re (C 



npq J 



[Ep — En){Eq — En) 



Here means that the terms with Ep — E^ are excluded. The rest are the terms where either p — n or q — n, and 



since Cnnq — ^nqn 



M\gl V\„, they contribute 

Cnnq lim (/„;„g(w, + fn-gn{^,^)) = ^ 



-3C„ 



{En - Eg)^ ■ 



Bringing together all the terms, we finally obtain the low-frequency limit for the quadratic susceptibility: 



Xo 



n _ n \ ^ M nq\ if^ nn qq) „ \ ^ \ ^ 

'-^^ {Eg-En)^ 



Re {C„ 



p q<p 



{Ep — En){Eg — En) 



XoA + XOB ■ 



(AlO) 



The meaning of the two terms in Eq. ljAlOp is different. Indeed, consider a system of noninteracting qubits. Then 



the eigenstates are completely factorized, \p) = 0^=1 \Pj)i ^^'^ the state of the jth qubit is either \gj) (ground) or \ej) 
(excited), with energy eg[e),j- Denoting by n the switched state {g — e,e ~ g), we obtain for Xqa 



M 

XoA,fact — 

i=l 



3\{n,\af\n^)f {{n,\af\ni) ~ {n,\(T^\n,)) 
(f- - f V ■ 



(All) 
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that is, the sum of single-qubit susceptibihties. On the other hand, one can check that 



is exactly zero for the system of independent qubits. Moreover, the summation restrictions in (jA10[l mean that the 
coefficients Cnpq can be replaced by the irreducible correlators (jlip of the operators erf: 

Cnp, M^C^^^j,^ - ^3"(^ J: k; n,p, q). (A13) 

Therefore the function Xob / (M^) will disappear in the thermodynamic limit in the absence of entangled 3-qubit 
clusters. Thus it provides an observable signature of 3-tanglement[23]. (Note that for = 3 the contribution 
from the disjoint correlators is exactly zero, due to cyclic invariance.) This signature can be extracted from the 
measurements by detecting the second harmonic generation when one, two etc qubits are near the co-degeneracy 
point and comparing the outcomes, much like it was done with the IMT deficit of |18| . 

Now we can introduce the iV-th order signature, Ea. (fT3l) . The functions Zjv are related to directly observable global 
nonlinear response of the system: the measurable quantity (xob)[13| is thus proportional to the same way the 
observed IMT deficit of Ref. [3 is proportional to Z2. We saw that the Zjy-functions are asymptotically zero in the 
thermodynamic limit, unless the system contains a finite proportion of entangled clusters of at least N qubits each; 
they are also directly related to the irreducible contribution to the corresponding static nonlinear susceptibilities. 
Therefore the functions (|13|14p can indeed serve as observable entanglement signatures. 
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